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We propose that cycle expansions be ordered with respect 
to stability rather than orbit length for many chaotic sys- 
tems, particularly those exhibiting crises. This is illustrated 
with the strong field Lorentz gas, where we obtain significant 
improvements over traditional approaches. 

PACS: 05.45.+b, 05.70.Ln 



Classical chaotic dynamical systems are inherently un- 
predictable, so that even if the inital conditions are 
known with high accuracy, predictions of the state of the 
system may be made to only a short time of order Aj~^ in 
the future, where Ai is the largest Lyapunov exponent. 
Despite this unpredictability, much can be said about the 
average behaviour of the system in many cases using only 
a small number of unstable periodic orbits or "cycles" . 
This can be achieved by using cycle expansions [0J|] of 
Ruelle's dynamical zeta function 

Cycle expansions have proved very useful in both clas- 
sical and quantum chaos, giving accurate estimates of the 
escape rate of open billiard systems and the energy 
levels of helium |^ using a surprisingly small number of 
classical cycles. The main idea of this approach is that a 
long generic trajectory may be approximated by various 
periodic orbits at different times, and that longer periodic 
orbits may often be "shadowed" by shorter "fundamen- 
tal" cycles, closely following the shorter cycles along dif- 
ferent sections of its length. Thus averages are calculated 
using fundamental cycles, with small corrections due to 
longer cycles. Periodic orbits which are exact repetitions 
of smaller cycles are explicitly summed, so that all ex- 
pressions are written in terms of the remaining "prime" 
cycles. The expansions work best when the the symbolic 
dynamics is well understood, and long periodic orbits are 
well shadowed by shorter ones. In this paper we investi- 
gate a system in which neither of these conditions holds, 
the strong field Lorentz gas. In spite of these difficul- 
ties, reasonable results may be obtained by ordering the 
expansion in terms of stability rather than the length of 
periodic orbits. This approach should be valid wherever 
cycle expansions can be applied, including flows for which 
a natural topological "length" is difficult to define. 

The origin of cycle expansions is a well developed the- 
ory of trace formulae and dynamical zeta functions . 
Here we will need only the expression for the classical 
time average of some quantity A in a closed system [|j , 
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The sum is over all distinct nonrepeating combinations 
of prime cycles, r is the period of a cycle, and A is mag- 
nitude of the expanding eigenvalue of the stability ma- 
trix, equal to exp(T E -^-i-)' where ^ A+ is the sum of the 
positive Lyapunov exponents. Note that for notational 
simplicity we have written A where other authors have 
|A|. Ai is the average of some quantity over a particular 
orbit, and could be Lyapunov exponent, current, or the 
number of collisions per unit time. An important check 
of the convergence of the expansion is given by 



C(0,0)-i = l + ^(-l)V(Ai---A,) 
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which should be equal to zero. 

The sums are usually truncated up to a particular 
order in the symbolic dynamics: for the Lorentz gas, 
the number of collisions N . At each order the terms 
are grouped into fundamental orbits that have a sym- 
bolic dynamics not separable into smaller cycles, and 
"curvature corrections" consisting of a composite or- 
bit, with symbolic dynamics ab together with its com- 
ponents a and b. Often the composite orbit has simi- 
lar averages to its components; then the combined term 
Kb + (AaTa + AbTb)A-^A^ ^ is smafi. If there 
are only a few fundamental orbits and this cancellation 
occurs most of the time the cycle expansion converges 
rapidly as a function of N. However, this is often not the 
case, as we see for the Lorentz gas (Tab. | below). 

We suggest that the orbit stability A may often be a 
more convenient parameter with which to truncate the 
cycle expansions. There are a number of reasons for 
this. The cycles which make the greatest contribution 
to the expansions are those with smallest A, which are 
not necessarily those of smallest N. The cycles of small- 
est A are also found most readily in numerical searches. 
Alternative methods for enumerating periodic orbits by 
exhaustive enumeration of symbol sequences (as in the 
zero field Lorentz gas |l^) fail whenever the symbolic 
dynamics is not understood completely and cycles exist 
which are not contained in the chaotic set. This is the 
case on one side of a chaotic crisis. See chapter 8 of 
Ref. [prl , where a number of examples are discussed, in- 
cluding the windows of the logistic map, and the Lorcnz 
attractor. Another reason for truncating in A rather than 
N is that there may be significant cycles (small A) at 
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comparatively large A''. Continuing the expansion to in- 
clude all cycles at large N may be prohibitive due to the 
large number of cycles, or that some of the cycles may 
have very large A and hence pose numerical problems. 
These issues are illustrated by our study of the strong 
field Lorentz gas. 

We now make the above prescription more precise. To 
retain the exponential convergence of traditional cycle 
expansions, it is clear that we must still make use of the 
close cancellation of curvature correction terms. This 
comes naturally if we include all terms for which the 
product of A's is less than a predetermined cutoff Amax- 
Then, if the product of A's match, the cancelling terms 
will be either included or excluded together. If the shad- 
owing is poor, cancellation would not occur whichever 
terms are included. In any case, it is clear that the ex- 
pansions will not work well if cycles are missing, whether 
the expansion is ordered with respect to N or A. Our 
prescription is thus: 

Alternative truncation of cycle expansions: only those 
cycles and combination of cycles which have the product 
of A 's less than a cutoff Amax should be included in the 
cycle expansion. Amax is the largest value such that al- 
most all the cycles are known. 

For the remainder of this paper we show how this works 
for the strong field Lorentz gas, obtaining significantly 
improved convergence. 

The two dimensional noncquilibrium Lorentz gas is a 
classical model of current flow in a conductor, and shares 
many of the properties of larger steady state noncqui- 
librium systems, including a simple relationship between 
the steady state current, and the sum of the Lyapunov 
exponents J • E = |l^JT^. Here, E is the imposed 
electric field. A single particle of charge e is scattered by 
a triangular lattice of circular disks of radius 1 and near- 
est neighbour spacing 0.236. Between collisions it moves 
according to the equations 



X = V V = eE/m — av 



(3) 



where a — ev-E/(mw^) is a Gaussian thcrmostatting co- 
efficient which keeps the kinetic energy constant, per- 
mitting the system to reach a steady current J — (ev)t. 
The symbolic dynamics is determined by the sequence of 
disks with which the particle collides; see Fig. |l|. 
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FIG. 1. The symbolic dynamics for the Lorentz gas at high 
field, for which only nearest neighbour collisions occur. We 
reduce the symbol set still further by noting that 10 is equiv- 
alent to 2 by symmetry, and similarly 8 is equivalent to 4. 



For small fields e = |eE/m| <C 1 the dynamics is known 
to be ergodic and lead to a steady current |[l5|]. There 
have been several previous investigations ]lO| , |l6| , |l7[ of the 
Lorentz gas at zero field (equilibrium) and small fields us- 
ing both cycle expansions (above) and direct calculation 
of measures based on periodic orbits. They used lO'^-lO^ 
cycles to evaluate the largest Lyapunov exponent, the 
diffusion coefficient (for zero field), and the current (for 
nonzero field). In each case the expansions were evalu- 
ated to a particular orbit length N, from 7 to 10. Most 
of the results have uncertainties of 5-10%, the main ex- 
ception being the evaluation of the Lyapunov exponent 
in Ref. |10|, which has an uncertainty of about 0.2%. 
This last calculation was performed in the fundamental 
domain, making full use of the symmetry of the lattice. 
The convergence for the Lorentz gas calculations is not 
as impressive as that in Refs. [^,^ because the grammar 
of the symbolic dynamics, that is, the rules determining 
which symbol sequences occur, is not well understood. 
Thus the amount to which longer orbits are shadowed by 
shorter ones is limited. 

At a field of e = 2.2 the ergodic attracting mea- 
sure is abruptly replaced by an attractor with fractal 
support characterised by a very restricted symbolic dy- 
namics |18 . In particular (see Fig. the only sym- 
bols are {2,4,6,8,10}, the allowed pairs of symbols are 
{26,28,48,410,62,64,68,610,82,84,104,106}, and the 
triples {262,264,1068,10610} are prohibited. These 
rules determine that the particle moves in one direc- 
tion along a channel formed by the disks, as described 
in Ref. [|l8|. Making use of the fact that the system 
is invariant under a reflection in the x-axis, we note 
that the symbols 2 and 10 are equivalent, and likewise 
4 and 8. Thus the above permitted sequences become 
{24, 26, 42, 44, 62, 64}, with no restrictions on the triples. 
It is easy to check that each symbol sequence of the 
restricted set corresponds to exactly two equivalent se- 
quences with the original symbols. Note that there are 
many other prohibited sequences of three or more sym- 
bols, depending on the field, so this system suffers from 
pruning in a similar fashion to the small field case. The 
range of fields we consider is 2.2 to 2.4, above which sta- 
ble orbits appear, rendering the attractor nonhyperbolic. 

The symbolic dynamics of this system is so restricted 
that, at e = 2.4 there are only nine orbits with seven 
or fewer collisions: see Tab. |. There are a couple of 
points to note from this table. First, there is very little 
shadowing of orbits. Five of the cycles are fundamen- 
tal in that they cannot be constructed from smaller cy- 
cles. In addition, the cycles that can be constructed from 
smaller cycles do not have eigenvalues A which are close 
to the product of the constituent A's, so curvature cor- 
rections obtained with these orbits are not particularly 
small. Secondly, the most significant orbits (with small 
A) do not occur at small A^. There are two = 10 orbits 
more significant than the N = 2 orbit, and 20 orbits with 
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A < 100, including one with = 34. Most of the longer 
orbits are partly cancelled in curvature corrections, but 
it is clear from the numerical results (below) that they 
are necessary to obtain optimal convergence. 

It is clear that obtaining all periodic orbits with fewer 
than, say, 30 collisions is quite unfeasible, since N — 7 
is already approaching the limit. In any case, this is 
not in line with the basic philosophy of describing the 
dynamics with as few cycles as possible. Thus we proceed 
to truncate the expansions using A rather than N. 

We numerically generated 10^ collisions for each of the 
field values from 2.2 to 2.4 in steps of 0.001, searching for 
periodic orbits of up to = 40. This generates about 
170 cycles for each field value, substantially fewer than 
Refs. ^0 16 1^], which considered a much smaller number 
of field and disk spacing values. Seven of the orbits in 
Tab. I were found, those with the smallest A. To get 
a more precise estimate of the A value at which cycles 
are first being missed, we plot the logarithm of the total 
number of cycles found as a function of In A in Fig. ||. The 
graph turns sharply down at about A = 600, indicating 
that the optimum A^ax should be around this value. 
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FIG. 2. The distribution of cycles found using our numer- 
ical routine, using bins of equal size in terms of In A. The 
maximum occurs at about A = 600, suggesting this as the 
optimal value of Amax- 

We computed C(0, 0)"^ using Eq. (Q) at each field value 
summing terms up to different values of N and Amax, and 
also evaluated the Lyapunov exponents Ai, A2 and the 
current J using Eq. (|l|) under the same conditions. The 
RMS averaged C(0,0)~^ values are shown in Fig. ||. The 
cycle expansion values for Ai, A2 and J were then com- 
pared with direct simulation results at each field value, 
to give the RMS differences exhibited in Tab. ||. 

From Fig. ^ we note that for each value of the cutoff, 
there is an initial regime in which convergence is exponen- 
tial with as usual for cycle expansions in hyperbolic 
systems [|[. At Amax = 600 this holds to about N = 24, 
while at Amax — 00 (no cutoff), convergence is expo- 
nential only to iV = 15. Subsequently the behaviour is 
determined by the limitations of the data set. If almost 
all of the cycles are present (as for A^ax < 600), the 



normalisation improves slightly, and then remains con- 
stant, otherwise it gets worse. A cutoff which is too 
high (Amax = 10'') can be more harmful than no cut- 
off, presumably because the unbalanced corrections at 
600 < A < 10^ are partly cancelled by the more numer- 
ous terms with opposite signs at larger A. 
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FIG. 3. Superior normalization of cycle expansions with 
stability cutoff, using Eq. . All expansions with Amax > 10 
converge exponentially to A*' = 15, after which the expansion 
with no cutoff (Amax ~ oo) diverges slightly, while the optimal 
cutoff (Amax ~ 600) continucs to improve up to A'^ = 30. 

The results for the Lyapunov exponents and the cur- 
rent (Tab. and Fig. m are more or less what would 
be predicted from the normalisation, except that the 
Amax = oo results are unexpectedly good, and compara- 
ble to Amax = 100. Given that J « 0.7, the relative error 
in the cycle expansion expression for J at Amax — 600 is 
about 0.3%, comparable to the best results in Ref. 
which required an order of magnitude more cycles. 




FIG. 4. The current at high field (center). Note the high 
degree of structure, which is reproduced in the cycle expan- 
sions with Amax ~ 600 (top) and to a lesser degree Amax ~ oo 
(bottom). The cycle expansion results are shifted by 0.05 for 
clarity. 

Finally, we use the results for ^(0,0)^' to estimate 
the rate of convergence of the expansion when almost all 
the cycles are known, that is, Amax < 600. Note that 
because a hyperbolic system has A ~ e^, exponential 
convergence in A'^ is equivalent to exponential conver- 
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gence in InAmax, that is a power law in Amax- This is 
shown in Fig. ||, where we find C(0,0)^^g ~ A^^^, or, 
C(0, 0)J^Jjg ^ 2"'"'^"'"=. The rate of convergence is thus 
comparable with that of cycle expansions performed up 
to a given N. The choice of which approach to take de- 
pends on whether it is more convenient to enumerate all 
orbits to a fixed iV or a fixed A, and the distribution of 

periodic orbits for the given system. 
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FIG. 5. The normalization at A'' = 40 as a function of 
In Aniax showing exponential convergence up to In Amax ~ 6.3. 
The solid line has a gradient of —0.69. 



It would be interesting to compare this method with 
the traditional approach for other classical systems, and 
also quantum systems, where formulae such as Eq. are 
replaced by more involved expressions, but still written 
in terms of classical periodic orbits. It could also be 
applied to flows for which no natural Poincarc section 
is known. This work was supported by the Australian 
Research Council. 
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r 


2 


26 


13.92 


1.459 


3 


264 


2.557 


2.706 


4 


2444 


5.337 


2.786 


5 


26264 


60.50 


4.133 


5 


24426 


175.6 


3.061 


6 


244426 


736.6 


4.147 


7 


2442644 


1.536 


5.763 


7 


2444264 


37.16 


5.471 


7 


2626264 


907.9 


5.594 



TABLE I. The shortest cycles at a field e = 2.4. Note the 
lack of shadowing, in that A244426 3> A2444A26. Also, A is 
not strongly correlated with N, so that truncating a cycle 
expansion at fixed A'^ tends to omit significant orbits with 
small A. 



-^max 


#av 


C(0,0)j^^s 


AAi RMS 


AA2 RMS 


A Jrms 


10 


3.7 


0.2283 


0.0559 


0.0865 


0.0249 


30 


10.2 


0.0647 


0.0648 


0.0774 


0.0136 


100 


26.8 


0.0292 


0.0226 


0.0312 


0.0061 


300 


62.9 


0.0153 


0.0095 


0.0125 


0.0029 


600 


104.0 


0.0134 


0.0074 


0.0102 


0.0025 


1000 


131.1 


0.0519 


0.0162 


0.0205 


0.0039 


10000 


157.6 


0.1640 


0.1814 


0.1525 


0.0256 


00 


169.2 


0.1264 


0.0190 


0.0331 


0.0083 



TABLE II. Cycle expansion results, showing the average 
number of periodic orbits used, #av, and the RMS differences 
of the other quantities from the direct simulation results. All 
contributions up to A' = 40 are included. The best conver- 
gence is obtained with Amax = 600, a factor of 3 better than 

no cut off (Amax = Oo) . 



